Spin-fluctuation dominated electrical transport of NisAl at high pressure 
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We present the first study of a magnetic quantum phase transition in the itinerant-electron fer- 
romagnet NiaAl at high pressures. Electrical resistivity measurements in a diamond anvil cell at 
hydrostatic pressures up to 100 kbar and temperatures as low as 50 niK indicate that the Curie 
temperature collapses towards absolute zero at a critical pressure pc — 82 ±2 kbar. Over wide ranges 
in pressure and temperature, both in the ferromagnetic and paramagnetic states, the temperature 
variation of the resistivity is found to deviate from the conventional Fermi-liquid form. We consider 
the extent to which this deviation can be understood in terms of a mean-field model of enhanced 
spin fluctuations on the border of ferromagnetism in three dimensions. 
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I. INTRODUCTION 

The electronic properties of metals on the border of 
magnetic phase transitions at low temperatures are of- 
ten found to exhibit temperature dependences at variance 
with the predictions of the standard model of a normal 
Fermi liquid. Early attempts to explain such non-Fermi- 
liquid behaviour have been based on a mean-field treat- 
ment of the effects of strongly-enhanced spin-fluctuations 
(paramagnons) iSi^ 

For a metal on the border of ferromagnetism in three 
spatial dimensions (3D), this mean-field spin- fluctuation 
model predicts a T^/^ temperature dependence of the 
resistivity^ instead of the conventional temperature 
dependence of a normal metal at low temperatures. The 
2^5/3 variation of the resistivity is a consequence of an 
underlying quasiparticle scattering rate that varies lin- 
early with the excitation energy of a quasiparticle near 
the Fermi level. This is the behaviour associated not 
with a Fermi liquid, for which the quasiparticle scatter- 
ing rate varies as E'^, but of a cross-over state known as 
the marginal Fermi liquid A review of several un- 

derlying models that yield a marginal Fermi-liquid form 
of the electron self energy is given in Ref. 0. 

In contrast to a phenomenological model introduced to 
describe the normal state of the cuprates;^ the marginal 
Fermi liquid that arises on the border of an itinerant- 
electron ferromagnet in 3D is due to the scattering of 
long-wavelength spin fluctuations that are relatively in- 
effective in reducing the current. This leads to a tempera- 
ture dependence of the transport relaxation rate or resis- 
tivity that is characterized by an exponent above unity, 
i.e., 5/3, but still below that expected for a conventional 
Fermi liquid^ 

In this paper we re-examine the behaviour of the resis- 
tivity on the border of itinerant-electron ferromagnetism 



in the relatively simple case of NisAl in which the cross- 
over from ferromagnetism to paramagnetism is achieved 
by the application of hydrostatic pressure. NiaAl can 
be prepared in a pure stoichiometric form and crystal- 
lizes in a simple cubic (CusAu) structured At ambient 
pressure it orders ferromagnetically below 42 K with a 
small average moment of 0.075 ^s/Ni in the limit of low 
temperature and low magnetic field>i2*ii*iSii^ 

For a test of the predictions of the mean-field spin- 
fluctuation model; NiaAl appears to offer advantages 
over a number of other metals on the border of mag- 
netism. The magnetic critical point can be tuned via 
hydrostatic pressure as opposed to chemical doping (as 
in, e.g., Pd(Ni))i^ that can introduce new physics of 
disorder not incorporated fully in the mean-field spin- 
fluctuation model. In contrast to other d-metals (e.g., 
MnSi, CoSi2)r^'^^'^^^ NiaAl has a more nearly contin- 
uous quantum phase transition and (unlike MnSi) docs 
not exhibit a spin-spiral structure due to the Moriya- 
Dsyaloshinski interaction arising from lack of space in- 
version symmetry. Also, in contrast to the nearly mag- 
netic f-metals (e.g., CcCug-xAuj;, CePd2Si2, YbRh2Si2, 
UGe2)r2iS242ii2S the spin-orbit interaction in NiaAl is rel- 
atively weak, the energy bands are relatively broad and 
the spin fluctuations are not local in real space, features 
that may be necessary for the applicability of the mean- 
field spin-fluctuation model in its present form. 

As a reference system for the study of the bor- 
der of metallic ferromagnetism, NiaAl is also con- 
venient because, along with its close relative, the 
nearly ferromagnetic metal NiaGa, it has been ex- 
tensively studied using a wide range of experimen- 
tal techniqueS,^-'^-^4i^5,26,27,28,29,30,31,32 ^j^j^j^ 

dicate the importance of spin fluctuations in these 
materialsi^i^ The combined results in these two systems 
define the parameters of the mean-field spin-fluctuation 
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model that may be relevant to interpreting the behaviour 
of Ni3Al from its low-pressure ferromagnetic phase to 
the high pressure paramagnetic phase near and above 
Pc where a marginal Fermi liquid cross-over state may be 
expected to be observedi^ 



II. EXPERIMENTAL DETAILS 

NisAl samples were pressurised using a diamond anvil 
cell (DAC) described in detail elscwherei^ The cell em- 
ployed anvils with culet diameter of 1 mm, a stainless- 
steel gasket with central hole for the sample and several 
ruby chips, and Argon as the pressure medium. After 
compression, the sample space was typically 400 /im in 
diameter and 50 /im in thickness. The pressure was de- 
termined from the fluorescence of the tiny ruby chips 
in the sample spacei^ The noble gas pressure medium 
ensured that the applied pressure was very hydrostatic. 
The pressure inhomogeneity over the sample was only 3% 
of the applied pressure. 

The samples were prepared from the melt by radio 
frequency heating starting with zone refined Ni and Al 
having residual resistivity ratios of over 2000. High ho- 
mogeneity samples were obtained by suitable stirring of 
the high purity stoichiometric melt, followed by rapid 
quenching and then annealing for up to 6 days. Small 
single crystals were spark cut from the resulting ingot 
and then characterised by means of microprobe analysis, 
transmission electron microscopy, Laue x-ray diffraction, 
mass spectroscopy, magnetic hysteresis and resistivity ra- 
tio measurements. These studies did not reveal evidence 
of precipitates of other phases or of total metallic impu- 
rity levels in excess of 20 ppm. The specimens selected 
for pressure studies in a ^He system and in a dilution 
refrigerator had residual resistivity ratios of 27 and 29, 
respectively, and had been used previously for the study 
of the de Haas- van Alphcn (dHvA) effect Tiny samples 
of the size required to fit into the DAC were cut to an 
initial thickness of 100 fim by means of low-power spark 
erosion and thinned to the final thickness of 10 fim by 
chemical polishing. 

The resistivity was measured by a sensitive ac 4- 
terminal technique. Four 12 /xm gold wires were at- 
tached to the sample by a micro-spot-welding technique 
that gave low contact resistance and thus low excitation- 
current heating and high detection sensitivity. Damage 
to the sample was minimized by spot welding with very 
low power immediately after the surface of the sample 
had been cleaned by chemical polishing. The gold wires 
were passed between one diamond anvil and one side of 
the stainless steel gasket insulated with a layer of 1266 
Stycast epoxy mixed with a saturated concentration of 
AI2O3 powder i2£ The resistivity was measured in a "^He 
system from 1.5 to 40 K with an excitation current of 
1 mA, and in a dilution refrigerator from 50 mK to 6.5 
K with an excitation current of 100 mA. The dilution re- 
frigerator contained a low temperature transformer and 
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FIG. 1: The temperature-dependent part of the electrical re- 
sistivity Ap = p — po at various pressures. Measurements were 
carried out in a ''He system (a) and a dilution refrigerator (b). 
Upon pressure application po increased irreversibly by about 
1 pQ,cn\ and thus the intrinsic pressure dependence of po could 
not be inferred. The curves have been shifted vertically for 
clarity. 



both refrigerators were equipped with a variable force ap- 
plication mechanism that allowed the in-situ changes of 
the pressure in the DACi^ 



III. RESULTS 
A. Magnetic phase diagram of Ni^Al 

Our measurements of the temperature dependence of 
the resistivity of NisAl at different pressures are shown 
in Figure nj The temperature dependent part of the re- 
sistivity is defined by Ap = p — po, where po is the resid- 
ual resistivity inferred by a suitable extrapolation of p 
vs T to absolute zero of temperature. The small val- 
ues of Po, initially of the order of 1 pficm for our sam- 
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FIG. 2: Signature of the ferromagnetic transition of NigAI in 
the temperature dependence of the resistivity. The ferromag- 
netic transition at Tcurie indicated by an arrow shows up in 
a plot of Ap — Apfit where Apfu is a smooth second-order 
polynomial fit of the data over the entire experimental range. 
The same values of Tcurie are obtained within experimental 
error from plots of dp/dT vs T. The curves have been shifted 
vertically for clarity. 
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FIG. 3: The proposed magnetic temperature-pressure phase 
diagram of Ni3Al. The full circles represent the peaks in the 
resistivity data Ap — Apju in Figure |5| The square repre- 
sents the ferromagnetic transition Tcurie = 41 ± 1 K at am- 
bient pressures inferred in samples similar to those studied 
here by means of both resistivity and magnetic susceptibil- 
ity measurements. The diamond defines the critical pressure 
Pc = 82 ± 2 kbar inferred from the peak in the T'^ coefficient 
of the resistivity vs pressure plotted in the inset of Figure Eb- 



pies, were found to increase irreversibly by about 1 fiflcm 
during the pressure experiments, and thus the intrinsic 
pressure dependence of po could not be inferred. At low 
and intermediate pressures, the resistivity curves exhibit 
very weak anomalies in p vs T centred at a tempera- 
ture that coincides at ambient pressure with the Curie 
temperature Tcurie- It is natural to associate the pres- 
sure dependence of these anomalies with the pressure de- 



pendence of Tcurie- The anomalies can be made more 
evident by plotting the derivative dp/dT or by sub- 
tracting from p a smooth polynomial fit of the form 
^Pfit = Co + ciT -I- C2T^- The curves Ap — Apfu shown 
in Figure 121 clearly reveal a peak that collapses monoton- 
ically with increasing pressure. Plots of dp/dT and of 
Ap — (co + ciT) also lead to peak positions consistent 
with those in Figure |21 The weakening of the peaks with 
increasing pressure is consistent with the predictions of 
the mean-field spin-fluctuation model. We note that near 
Tcurie the magnetic contribution to dp/dT is quite gen- 
erally expected to be proportional to the magnetic con- 
tribution to the heat capacity^S*^ which, in the mean- 
field spin-fluctuation model, is predicted to decrease with 
decreasing Tcurie^ Particularly at high pressures, how- 
ever, the peaks may also be reduced by broadening due 
to pressure inhomogeneities. This effect grows with in- 
creasing downward curvature of Tcurie vs p and may be 
expected to lead to a disappearance of the Curie point 
anomaly in p vs T near pc- 

The pressure dependence of Tcurie inferred from the 
resistivity anomalies is shown in Figure |21 Tcurie col- 
lapses with increasing pressure with an initial gradient of 
approximately — 0.4±0.05 Kkbar"^, which is close to the 
value reported in an early low-pressure studyi^ Assum- 
ing a slightly stronger than linear variation, Tcurie vs p 
extrapolates to a value close to Pc ~ 82 ± 2 kbar defined 
by the peak in the low-temperature limit oi A = Ap/T^ 
vs p (see Figure [71l. However, the lack of evidence of a 
clear signature of Tcurie in Ap vs T above 60 kbar leaves 
some doubt as to the true value of the critical pressure 
at which Tcurie vanishes. 

We note that ferromagnetism in NiaAl can also 
be suppressed by Al doping piii^ However, the high 
sensitivity of Tcurie to dopant concentration x in 
Ni75_2;Al25+2; has made it difficult to produce a detailed 
temperature-dopant phase diagram for comparison with 
our temperature-pressure phase diagram. In previous 
work X was varied in steps of 0.5, corresponding to steps 
of the order of 100 kbar, and the critical value of x where 
Tcurie vauishcs is estimated by interpolation of the data 
to be about 0.4. Doped specimens have higher po and 
thus lower values of Ap/po, making measurements of the 
temperature dependence of p more difficult. Also, doped 
samples may have significant inhomogeneities in Tcurie 
and may involve physics not included in the mean-field 
spin-fluctuation model. 



B. Temperature dependence of the resistivity 

The evolution with pressure of the temperature de- 
pendence of the resistivity is shown in Figures ^ and |S1 
The resitivity is plotted against T, T^/^, T^/a^ j.2 
for three pressures corresponding to the lower and upper 
end of the pressure range of the data in Figure ^3 as well 
as the critical pressure. 

Overall the temperature dependent part of the resistiv- 



4 



ity appears to be only weakly pressure dependent over the 
entire pressure range investigated up to nearly 100 kbar. 
This is in marked contrast to the behaviour of MnSi and, 
in particular, of narrow f-band materials that in general 
exhibit obvious and strong variations of the resistivity 
upon crossing the critical pressure over a pressure range 
of typically only a few kbar. NisAl differs in its low- 
temperature resistivity from typical f-band systems in 
being essentially very near a magnetic quantum critical 
point over a wide pressure range. 

Figures El and also show that even at low tempera- 
tures, 0.05 K to 7 K, the resistivity cannot be described 
in terms of a simple power law. The overall best fit in 
this temperature range is to an exponent of the order 
of 3/2. Below 4 K, however, the best fit is to an expo- 
nent of the order of 5/3 and below 1 K to an exponent 
of about 2. This behaviour is highlighted in the tem- 
perature variation of the logarithmic derivative of the re- 
sistivity 91nAp/91nr shown in Figure El which defines 
the temperature dependent resistivity exponent n{T). At 
still higher temperatures (well above 10 K), the resistivity 
exponent eventually drops towards unity (see Figure^). 

For low pressures where comparisons can be made our 
findings are generally consistent with previous resistivity 
measurements in stoichiometric sampleaiS'iS'Sii^il^ and 
in samples in which the ratio of the Ni and Al concentra- 
tions were varied so as to suppress ferromagnetic order^ 
In all cases the temperature dependence of the resistiv- 
ity is found to deviate from the Fermi liquid form and 
to be described in first approximation in terms of a re- 
sistivity exponent around 3/2 in the temperature range 
1-30 K. However, in contrast to the present work, these 
earlier studies did not yield clear evidence for a limit- 
ing form of the resistivity at very low temperatures 
below 1 K. In particular, Fluitman et al., found a resis- 
tivity exponent of the order of 3/2 down to 200 mK in 
samples with Ni to Al ratios tuned to the critical value 
where Tcurie vanishes 4& Furthermore, in a study of sev- 
eral stoichiometric NisAl samples, Stciner et al., found 
that while above 1 K the resistivity exponent is sample 
independent, and consistent with our work, below 1 K 
it varied significantly from sample to sample, exhibiting 
in some cases weak upturns or downturns below a few 
hundred mKi^ 

We note that in all samples investigated thus far Ap 
is very small compared with po at dilution refrigerator 
temperatures (i.e., Ap/po < 2% for T < 1 K). Thus, 
the behaviour of Ap in this range may be sensitive to 
sample heterogeneities and the precise way by which the 
samples are prepared. Thus, the form of the resistivity 
that we observe below 1 K (Figure [TJ is not necessarily 
the property of the ideal stoichiometric state of NisAl. 
It is interesting to note, however, that the coefficient 
A of Ap in the mK range shows a peak (the peak po- 
sition we define to be Pc, inset of Figure E^i). Plots of 
Ap/T^ vs T at different pressures in Figure |H1 show that 
the effect of pressure and the approach of a quantum crit- 
ical point near pc become most clearly evident in the low 
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FIG. 4: Temperature-dependent part of the resistivity Ap 
of NisAl versus T (a) and versus T^^^ (b). Ap varies more 
strongly than linear with temperature. At low temperatures, 
Ap rises more strongly than T^^^ and at higher temperatures 
Ap rises more weakly than T"^/^. As shown in Figure El the 
resistivity exponent, defined by the logarithmic derivative of 
Ap, is not constant and varies strongly with temperature even 
in the liquid helium temperature range. 



temperature limit. We see that not only is the maximum 
of Ap/T^ vs T at Pc strongest at low temperatures but 
also the crossover temperature above which p clearly de- 
viates from the Fermi- liquid form (e.g., the temperature 
where Ap/T^ falls by 20% of its value in the zero tem- 
perature limit) reaches a minimum at this same pressure. 
At 81 kbar Ap/T^ vs T continues to grow with decreas- 
ing temperature down to at least below 600 mK where 
^ Po and a well-defined finite value of ^ = Ap/T^ 
as given in the inset of Figure \7}p is obtained only after 
averaging over 500 mK. Although Fluitman et al. did 
not report observing a form of Ap, they nevertheless 
considered the behaviour of the average value between 
1.2 and 4.2 K of the ratio Ap/T'^M They find that this 
ratio varied by a factor of about 1.3 in going from a sto- 
ichiometric NisAl sample to a non-stoichiometric sample 
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FIG. 5: Temperature-dependent part of the resistivity Ap 
of NisAl versus T^^^ (a) and versus (b). The resistiv- 
ity exponent defined by the logarithmic derivative of Ap is 
not constant and increases with decreasing temperature from 
about 1.5 at 5 K towards 2 below 1 K (Figure El- 



with critical doping where Tcm-ie vanishes. This agrees 
well with the variation versus pressure of ^pjT'^ in the 
same temperature range (Figure IHJ. 



IV. DISCUSSION 

In the following we discuss our experimental results 
in terms of an elementary form of the mean-field spin- 
fluctuation model discussed in the introduction^^ Near 
the critical pressure where ferromagnetic order in 3D van- 
ishes, this model leads to an electron self energy of a form 
at low energy characteristic of a marginal Fermi liquid 4^ 

We consider a 3D isotropic itinerant-electron system 
in which the total spin is conserved and only long wave- 
length fluctuations of the magnetization are strongly en- 
hanced by an exchange molecular field arising from the 
effects of the Coulomb interaction and the Pauli princi- 
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FIG. 6: The resistivity exponent of NisAI defined by the log- 
arithmic derivative of the resistivity, n = dlnAp/ dlnT . The 
resistivity exponent rises with decreasing temperature from 
about 1.5 at 5 K towards 2 below 1 K. At each T the ex- 
ponent is a minimum at around p^. We note that there in 
no extended temperature regime for either a or T^^^ be- 
haviour. The error bar given at 5 K applies to all data points 
shown in the plot. 



pie. These conditions may be approximately satisfied in 
the case of NiaAl for the reasons given in the introduc- 
tion. The essential features of the model are as follows. 
We assume that in the low T limit the magnetization 
M(r) in a weak applied magnetic field H(r) is given by 
a Ginzburg-Landau equation of the form 



H = aM + hm^ - cV^M 



(1) 



where M'^ — (M • M)M and h and c are positive con- 
stants. The latter conditions imply that the system un- 
dergoes a continuous ferromagnetic transition when a 
crosses zero. As in the Landau mean-field model we as- 
sume a is linear in {jp — pc) for p close to the critical 
pressure pc, i-e., that a = a(j) — pc), where a is a positive 
constant. An analytic expansion of a and the mean-field 
analysis given below is thought to be plausible because 
the effective dimension relevant to quantum phenomena 
in the T —> limit is greater than the upper critical di- 
mension in NiaAliSi^ We also assume that relaxation of 
a fluctuation of the magnetization M(r, t) to equilibrium 
is governed by Landau damping. Thus, a Fourier compo- 
nent Mq(i) of M(r,t) decays exponentially to the value 
given by Equationnvia a relaxation function that is pro- 
portional to q = |q|. In zero applied field and small Mq 
in the paramagnetic state this implies 
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-jq{a + cq^)Mq 



(2) 



where 7 is a positive constant. We note that the coef- 
ficient of (^Mq) on the right hand side of Equation El 
defines the q dependent relaxation rate Tq of magnetic 
fluctuations. Thus, Tq is given by 7(7 times the inverse 
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FIG. 7: Temperature dependent part of the resistivity Ap of 
NisAl in the millikelvin range. In this range where Ap is only 
of the order of 1-2% of po, Ap is approximately of the form 
AT'^ for our sample. The coefficient A increases with rising 
pressure up to pc = 82 ± 2 kbar (a) and falls again beyond pc 
(b). As shown in the inset of Figure (b) the peaks value of 
A at Pc is a factor of 2 higher than the zero-pressure value of 
0.01 pf2cmK^'^^ A has been determined from a fit of Ap vs 
in the range 50 mK to 600 mK. 



static susceptibility, which is defined as the linear coef- 
ficient of Mq in the Fourier transform of Equation ^ 
Note that when a — > the relaxation rate becomes cubic 
in q, i.e., the dynamical exponent z is here equal to 3. 
Thus for 13 = 3 the effective dimension D ~\- z relevant 
to quantum critical phenomena is equal to 6, which is 
indeed above the upper critical dimension of 4 for our 
Ginzburg-Landau model. 

The model defined by Equations ^ and |21 is found 
to be consistent with ambient-pressure bulk magne- 
tization measurementsii that yield the parameters a 
and h and ambient-pressure inelastic neutron scattering 
measuremcntai^ that yield estimates of c and 7 in NisAl. 
The values of these essentially ground state parameters 
define the starting point of a model for the temperature 
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FIG. 8: Pressure dependence of Ap/T^ versus T in NiaAl. 
The ratio /^p/T"^ increases with decreasing temperature and 
tends to saturate at low temperature. The characteristic tem- 
perature at which saturation occurs decreases with increasing 
pressure up to the critical pressure (p^ = 82 ± 2 kbar). Cor- 
respondingly, the saturation value of /^p/T"^ increases with 
pressure and peaks near pc (see inset of Figure (TJa). We note 
that near pc the saturation of IS.pjT'^ is not evident except 
when averaging over the data from 50 to 600 mK (squares) 
where Ap <C po. 



dependence of the magnetic equation of state, heat ca- 
pacity and resistivity. 

The temperature dependence of the magnetic equation 
of state is assumed to arise primarily from the effects 
of strongly enhanced long-wavelength spin fluctuations 
rather than from incoherent particle-hole excitations of 
conventional band theory. Specifically, any temperature 
dependence is assumed to arise in a self-consistent way 
from Equations n and 121 together with the Bose function 
which governs the excitations of those spin fluctua- 
tions that arc strongly excited on the border of magnetic 
long-range order. Thus, the finite temperature proper- 
ties are consequences of a set of T = parameters and a 
universal thermal factor in analogy to the Debye model 
for lattice vibrations or the Fermi liquid model for the 
quasiparticle excitations of a normal metal at low tem- 
peratures. The principal difference is that in the spin- 
fluctuation model the relevant modes that are thermally 
excited are not normal modes in the usual sense, but are 
characterized by a relaxation spectrum (see Equation|2Jl. 
The self-consistent spin-fluctuation model and the ear- 
lier non-self-consistent paramagnon model can thus be 
viewed as being elementary examples involving the sta- 
tistical mechanics of open or dissipative systemsi^2iSi*^ 



A. Magnetic phase boundary 

We consider flrst the effect of thermal fluctuations 
of the magnetization on the magnetic equation of state 
(Equation ^ . These fluctuations may be imagined to 
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arise from the effect of a random field of zero mean added 
to the left hand side of Equation ^ which then drives a 
random magnetization of zero mean added to the aver- 
age magnetization in each term on the right hand side 
of Equation ^ In the absence of the non- linear term on 
the right hand side of Equation^such fluctuations when 
averaged over an ensemble yield no effect. However, the 
non-linear term (feM'') leads to corrections that depend 
on the variance of the local magnetization. In lowest 
order, in the paramagnetic state, for example, one finds 
that the linear coefhcient in Equation^ which represents 
the inverse uniform susceptibility, becomes 



m 



(3) 



(4) 



(5) 



where ^|mq| J is the thermal variance of a Fourier com- 
ponent of the fluctuating component of the magnetiza- 
tion and the sum is per unit of volume. The thermal 

variance ^|mq|^^ is defined by the fluctuation-dissipation 

theorem in terms of ji^^ and the wavevector and frequency 
dependent susceptibility Xqw In our model the latter is 
given by Equations Q and 13 but with a replaced as in 
Equation This yields a self-consistent equation for 

^|mq|^^ and thus the magnetic equation of state versus 

temperature and pressure. (Note that the zero-point con- 
tribution to the total variance of the local magnetization 
is not included in Equation [3 For further discussions 
see, e.g., Refs. 31 and 43.) 

The Curie temperature is defined by a vanishing lin- 
ear coefficient in the magnetic equation of state. Equa- 
tions then yield under our assumption and for p < pc, 



2.39C73 



a(Pc - p) 



(6) 



Therefore, near pc, the Curie temperature should vary 
as {pc — pY^'^ and depend solely on the ground state 
parameters of the model defined by Equations ^ and |21 
This result holds only for sufficiently low values of Tcurie 
where the breakdown region of the mean-field approxi- 
mation around Tcurie is sufficiently narrow. From the 
measured parameters Equation yields a value of Tcurie 
of the order of 40 K for NiaAl at ambient pressure, in 
good agreement with experiment Our observed pres- 
sure dependence of Tcurie is not inconsistent with the 
model, although more detailed measurements near pc are 
needed to provide a test of the predicted ipc—pY^'^ vari- 
ation of Tcurie- The phcnomcnological parameters in 



the model might be inferred from an appropriate energy- 
band model that includes effects of zero-point spin fluctu- 
ations. Recently, such a model has been considered and it 
leads in particular to insights into the origin of the pres- 
sure dependence of the ordered moment in NiaAli^iM 



B. Temperature dependence of the resistivity 

Next we consider the effect of thermal spin fluctuations 
on the resistivity. Within the Boltzmann-Born modelS 
one finds that Ap in the paramagnetic state can be ex- 
pressed in the form 




(7) 



where 77 is a constant. The temperature derivative is 
solely with respect to the Bosc function entering the def- 
inition of the thermal variance. We have assumed that 
momentum is efficiently transferred from spin fluctuation 
to the lattice either by Umklapp processes or residual dis- 
order. The thermal variance of the local magnetisation 
arises in Equation[71because current carriers are assumed 
to undergo weak scattering from a molecular exchange 
field that is proportional to the local magnetization. Mul- 
tiple scattering processes are not included. The factor q 
in Equation[71comes from a product of (i) a factor aris- 
ing from the fact that low q fluctuations are ineffective 
in reducing the current and (ii) a partly compensating 
\/ q factor due to a loss of wavevector phase space com- 
ing from the Pauli principle that constrains scattering 
to the vicinity of the Fermi surface. Corrections due to 
momentum non-conservation possibly arising from resid- 
ual disorder or intraband transitions are not included. 
The temperature derivative Td/dT can be shown to arise 
essentially from the fact that the available momentum 
phase space increases with increasing temperature due 
to the thermal depopulation of states below the Fermi 
surface. 

Below a characteristic temperature Tfl the above 
model (Eauations lll7|l leads to Ap of the form Ap = AT^ 
where A is inversely related to the square of the tempera- 
ture TpL that vanishes as p — > Pc- At this pressure Ap is 
predicted to vary as T^/^ in the low T limit. An elemen- 
tary discussion of (i) the connection of Tp l to the model 
parameters and (ii) the T^/^ form of the resistivity and 
its connection to the marginal Fermi liquid model may 
be found for example in Ref . ITfil 

An extension of the above treatment into the ferro- 
magnetic state^*^ leads again to a resistivity below 
a characteristic temperature TpL and a T^/^ resistivity 
above Tcurie, but below another characteristic temper- 
ature Tmfl- Above Tmfl the resistivity exponent de- 
creases and tends to unity at higher T due to the growth 
of the thermal variance of the local magnetization in 
Equation 131 i.e., due to the coupling between spin fluc- 
tuation modes in the mean field approximation. 
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FIG. 9: Temperature-pressure phase diagram predicted by 
the mean-field spin-fluctuation model discussed in the text. 
For p > Pc the lower axis represents the square of the mag- 
netic correlation vector (inverse magnetic correlation length) 
in the zero temperature limit. The solid line represents the 
Curie temperature given by Equation |S| The crossover lines 
Tfl and Tmfl are defined in the text. If the ferromagnetic 
transition becomes first order at low temperatures, then there 
exists a region (vertical shading) of forbidden values of a/ci— 



C. Temperature-pressure phase diagram for NiaAl 

The above findings are summarized in the predicted 
temperature-pressure phase diagram for NiaAl shown in 
Figure 13 The pressure dependences of Tcurie below Pc 
and of Tmfl and Tfl in the paramagnetic state below 
or above Pc are obtained from (i) Equations I1I7I (ii) 
the model parameters as defined in Ref. Issl and (iii) 
Pc = 82 kbar. The values of TpL below pc are ob- 
tained by an extension of the model to the ferromag- 
netic statei^*^ The ferromagnetic transition is assumed 
to be everywhere continuous. The Curie temperature 
Tcurie is defined by Equation |3 The characteristic 
crossover temperature TpL is defined by the condition 
n{TpL) = InAp/lnT = 1.8 and the crossover tempera- 
ture Tmfl by the condition n{TMFL) ~ 1-6. 
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FIG. 10: The temperature dependent part of the resistivity 
Ap vs T^l'-^ normalized to values at Tcurie- Predictions of the 
mean-field spin-fluctuation model (a) and experimental find- 
ings (b) for two samples with different values of Tcurie- The 
higher value of Tcurie was obtained by chemical dopingS^ and 
the lower value by the application of hydrostatic pressure to 
our stoichiometric samples (Figure^). The correspondence 
between theory and experiment is strikingly close except for 
the sample with the lower Tcurie at low temperatures. This 
discrepancy at low T, which indicates an extra component 
in Ap additional to the spin- fluctuation component Apsf, ex- 
tends up to Pc and beyond. It is also visible in earlier ambient 
pressure data in a stoichiometric sample^^*^^*^ 



D. Comparison with experiment 

We find that the experimental results may be under- 
stood quite well in terms of the above model over a wide 
range in temperature and pressure, except in the range 
of about 1-10 K. The regions of agreement and disagree- 
ment between the model and experiment are illustrated 
in the examples shown in Figure [TUl The variation of Ap 
vs T^/^ in two samples of NisAl with different Tcurie is 
plotted in normalized form in Figure [TUb ^ The predic- 
tions of the model for the set of parameters appropriate 
to the two samples are shown in FigurclTUki^ The agree- 
ment between theory and experiment is rather striking 
both in the magnitudes of Tcurie and the overall form of 



the curves. 

However, the stoichiometric sample at high pressure 
exhibits an anomalous downturn at low T which is not 
anticipated by the model. This discrepancy is present 
in all of our high precision measurements in stoichiomet- 
ric samples up to Pc and beyond. It is also evident in 
an earlier study at ambient pressured The above model 
predicts a T^/^ form of Ap at all T below approximately 
10 K at the critical pressure. This is not observed (Fig- 
ure (H)). Instead, as discussed earlier, we observe a con- 
tinuous change in x{T) varying from approximately 1.3 
towards 2 as the temperature is decreased from 5 K to 
below 1 K. 
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The temperature dependence seen away from pc is 
qualitatively consistent with the model. However, the 
magnitude of Tp^ estimated from Figure El or from Fig- 
ure IHl appears to be nearly an order of magnitude lower 
than predicted by the model (Figure ISJ. 



E. Corrections to the elementary mean-field 
spin-fluctuation model 

The existence of a regime at pc well below 1 K is not 
expected if the ferromagnetic transition is second order. 
A finite Tp l at pc may indicate that the phase transition 
becomes first order close to pc as might be generally ex- 
pected for ferromagnetic quantum phase transitionsi^S^. 
Here wc note that, as found in MnSi, a first order tran- 
sition need not lead to a Fermi liquid form of Ap near pc 
(i.e., where Tcurie vanishes) 

The most significant challenge to the mean-field spin- 
fluctuation model is to account for the anomalous down- 
turn seen in Figure ITnb . which was not anticipated and 
remains unexplained. We stress, however, that in other 
respects the mean-field spin-fluctuation model provides a 
rather accurate description of the temperature-pressure 
phase diagram of NisAl. 

The low temperature form of the resistivity is also 
poorly understood in a number of other materials. In 
MnSi one observes a T^l'^ variation of Ap over nearly 3 
decades below 10 K in an extended range in pressure 
beyond the critical pressure of a first order magnetic 
transitioni^SiffiiSi This is at odds with the prediction that 
Tp L should be of the order of 5 K or higher at all pres- 
sures in this system. In e-Fe the resistivity varies as T^/^ 
over two decades above the critical pressure of a first or- 
der magnetic and structural phase transition. Naively 
one would have expect to see a resistivity below a few 
K in this system too^^^ 

Non- Fermi liquid forms of the resistivity have now been 
observed in a large number of examples in which the 
applicability of the mean-f ield sp in-fluctuation model is 
questionable (see, e.g., Refs l45lMl65f) . The cases of MnSi 
and NiaAl may stand out partly because the magnitudes 
of the discrepancies between the predictions and the ex- 
perimental findings can be quantified in terms of inde- 
pendently determined model parameters. Also, in these 
two materials the discrepancies arise in small pockets of 
the temperature-pressure phase diagrams surrounded by 
wide regions where the mean-field spin-fluctuation model 
provides a rather accurate description of the experimen- 
tal findings. 

The systems in which the discrepancies between theory 
and experiment arc most clearly evident are character- 
ized by first order transitions for sufficiently low Tcurie- 
This implies that the parameter b in Equation ^ must 
be taken to be negative, i.e., the coupling between the 
spin-fluctuation modes is attractive. The applicability 
in this case of the mean-fleld decoupling method is not 
self-evident. Systems described by a quantum Ginzburg- 



Landau model in which & < (but higher order coeffi- 
cients in the expansion in M are positive) may be charac- 
terized by highly anharmonic magnetic fluctuations not 
described in the mean-field approximation. Such fluctu- 
ations may be of particular interest in cases when the 
tricritical point characteristic of this extended quantum 
Ginzburg-Landau model is at sufficiently low tempera- 
tures. 

Beyond this one may also look for an explanation 
of the discrepancies between experiment and the mean- 
field spin-fluctuation model in a number of more con- 
ventional effects arising from lattice vibrationSfi2iS& band 
degeneracies |Si Umklapp processeaS^ and quenched disor- 
der. A systematic study of such effects is still lacking for 
the systems we have considered. We note, however, that 
the effects of conventional phonon scattering, in particu- 
lar, is expected to be very weak and ignorablc in NiaAl 
at least below 100 KiiS Also, the effects of quenched dis- 
order may be expected to be strongest below 1 K where 
Ap <^ po in our NiaAl samples. However, the behaviour 
of Ap vs T in this regime is not inconsistent with our 
model predictions except perhaps very near pc- For ex- 
ample, Ap/T^ at 500 mK is predicted to increase by a 
factor of about 2 from ambient pressure to p^ ■ This agrees 
well with the variation of A shown in the inset of Fig- 
ure [7|3. These considerations suggest that the anomalous 
behaviour of Ap in the range 1-10 K may arise from essen- 
tially intrinsic effects of the electron-electron interaction 
and not explicitly from phonons or quenched disorder. 



V. CONCLUSION 

The magnetic temperature-pressure phase diagram for 
NiaAl has been explored for the first time beyond the 
low-pressure regime by means of 4 terminal resistivity 
measurements in a diamond anvil cell. Hydrostatic pres- 
sure studies up to 100 kbar and down to 50 mK indicate 
that the Curie temperature collapses towards absolute 
zero at a critical pressure Pc = 82 ± 2 kbar. 

A non-Fermi liquid form of the resistivity is observed 
over a wide range in temperature and pressure except be- 
low a characteristic temperature Tp l that decreases from 
a few Kelvin at ambient pressure to well below one Kelvin 
at Pc- The finite value of Tp^ at pc may indicate that 
the ferromagnetic transition is first order near the critical 
pressure. The temperature dependence of the resistivity 
and Tcurie is, overall, consistent with the predictions of 
the mean-field spin-fluctuation model, which reduces to 
a type of marginal Fermi liquid model just above Tcurie 
when Tcurie tends towards absolute zero. However, a 
comparison of theory and experiment also reveals the ex- 
istence of an extra component in the temperature depen- 
dence of the resistivity at intermediate temperatures that 
has a different and unidentified origin. 



10 



Acknowledgments is grateful for support by the FERLIN program of the 

European Science Foundation. 

We wish to thank S. Julian, A. Rosch, A. Huxley and 
S. Yates for stimulating discussions on this topic. PGN 



* e-mail: niklowit Qdrf mc . ceng . cea. t'r present address: 
Departement de Recherche Fondamentale sur la Matiere 
Condensee, SPSMS, CEA Grenoble, 38054 Grenoble 
Cedex 9, France 

^ J. A. Hertz, Phys.Rev.B 14, 1165 (1976), and references 
therein. 

^ A. J. Millis, Phys.Rev.B 48, 7183 (1993), and references 
therein. 

^ T. Moriya, Spin Fluctuations in Itinerant Electron Mag- 
netism (Springer, Berlin, 1985), and references therein. 

* J. Mathon, Proc.RoySoc.A 306, 355 (1968). 

^ C. M. Varma, P. B. Littlewood, S. Schmitt-Rink, E. Abra- 
hams, and A. E. Ruckenstein, Phys. Rev. Lett. 63, 1996 
(1989). 

^ T. Holstein, R. E. Norton, and P. Pincus, Phys.Rev.B 8, 
2647 (1973). 

G. Baym and C. Pethick, Landau-Fermi liquid theory (Wi- 
ley, New York, 1991), chap. 3. 

^ A relaxation rate linear in E might be expected to give rise 
to a resistivity varying as T rather than as T^/^ , The extra 
factor of r^/^ arises in our case because spin-fluctuations 
or paramagnons on the border of ferromagnetism give rise 
to scattering mainly in the forward direction and are thus 
relatively ineffective in reducing the current. 

® W. G. Moffatt, The Handbook of Binary Phase Diagrams 

(Genium Publishing Corporation, New York, 1983). 
^° J. H. Fluitman, B. R. de Vries, R. Boom, and C. J. 
Schinkel, Phys.Lett.A 28, 506 (1969). 

F. R. de Boer, C. J. Schinkel, J. Biesterbros, and S. Proost, 
J.Appl.Phys. 40, 1049 (1969). 
^"^ H. Sasakura, K. Suzuki, and Y. Masuda, J.Phys.Soc.Jpn. 
53, 352 (1984). 

N. R. Bernhoeft, Ph.D. thesis. University of Cambridge 
(1982). 

" M. Nicklas, M. Brando, G. Knebel, F. Mayr, W. Trinkl, 
and A. Loidl, Phys.Rev.Lett. 82, 4268 (1999). 

^5 C. Pfleiderer, G. J. McMullan, S. R. Juhan, and G. G. 
Lonzarich, Phys.Rev.B 55, 8330 (1997). 

C. Thessieu, C. Pfleiderer, A. N. Stepanov, and J. Flou- 
quet, J.Phys.Cond.Mat. 9, 6677 (1997). 

" N. V. Mushnikov and T. Goto, Phys. of Metals and Met- 
allography 93, S88 (2002). 

S. Barakat, Ph.D. thesis, University of Cambridge (2004). 
^® A. Schroder, G. Aepph, R. Coldea, M. Adams, O. Stockert, 

H. Lohneisen, E. Bucher, R. Ramazashvili, and P. Cole- 
man, Nature 407, 351 (2000). 

J. Custers, P. Gegenwart, H. Wilhelm, K. Neumaier, 
Y. Tokiwa, O. TrovareUi, C. Geibel, F. Steglich, C. Pepin, 
and P. Coleman, Nature 424, 524 (2003). 
S. S. Saxena, P. Agarwal, K. Ahilan, F. M. Grosche, 
R. K. W. Haselwimmer, M. J. Steiner, E. Pugh, I. R. 
Walker, S. R. Juhan, P. Monthoux, et al.. Nature 406, 
587 (2000). 

22 N. D. Mathur, F. M. Grosche, S. R. Julian, I. R. Walker, 

D. M. Freye, R. K. W. Haselwimmer, and G. G. Lonzarich, 



Nature 394, 39 (1998). 

J. H. J. Fluitman, R. Boom, P. F. de Chatel, C. J. Schinkel, 
J. L. L. Tilanus, and B. R. de Vries, J.Phys.F 3, 109 (1973). 
2'' N. Buis, J. J. Franse, and P. E. Brommer, Physica B 106, 
1 (1981). 

2^ T. I. Sigfusson, N. R. Bernhoeft, and G. G. Lonzarich, 
J.Phys.F:Met.Phys. 14, 2141 (1984). 

J. J. M. Buiting, J. Kuebler, and F. M. Mueller, J.Phys.F 

13, L179 (1983). 
2^ N. R. Bernhoeft, G. G. Lonzarich, P. W. MitcheU, and 

D. M. Paul, Phys.Rev.B 28, 422 (1983). 
2* N. R. Bernhoeft, G. G. Lonzarich, P. W. MitcheU, and 

D. M. Paul, PhysicaB 136, 443 (1986). 
28 N. R. Bernhoeft, S. M. Hayden, G. G. Lonzarich, D. M. 

Paul, and E. J. Lindley Phys.Rev.Lett. 62, 657 (1989). 
^° S. K. Dhar, K. A. Gschneidner, L. L. Miller, and D. C. 

Johnston, Phys.Rev.B 40, 11488 (1989). 

F. Semadeni, B. Roessli, P. Boni, P. Vorderwisch, and 
T. Chatterji, Phys.Rev.B 62, 1083 (2000). 

^2 S. H. Kilcoyne and R. Cywinski, PhysicaB 326, 577 (2003). 
A. Aguayo, L I. Mazin, and D. J. Singh, Phys.Rev.Lett. 
92, 147201 (2004). 

L L Mazin, D. J. Singh, and A. Aguayo, cond- 
mat/0401563. 

G. G. Lonzarich and L. Taillefer, J.Phys.C:Solid State 
Phys. 18, 4339 (1985). 

J. Thomasson, F. Thomas, C. Ayache, L L. Spain, and 
M. Villedieu, in Frontiers of High Pressure Research, 
edited by H. D. Hochheimer and R. D. Etters (Plenum, 
1991), pp. 423-432. 

A. Jayaraman, Rev.Mod.Phys. 55, 65 (1983). 
'"^^ J. Thomasson, Y. Okayama, I. Sheikin, J.-P. Brison, and 

D. Braithwaite, SolidStateComm. 106, 637 (1998). 
^8 B. Sake, J. Thomasson, A. Demuer, J. J. Blanchard, J. M. 

Martinod, L. Devoille, and A. Guillaume, Rev. Sci. lustrum. 

71, 2461 (2000). 
""^ M. E. Fisher and J. S. Langer, Phys.Rev.Lett. 20, 665 

(1968). 

S. Alexander, J. S. Helman, and I. Balberg, Phys.Rev.B 
13, 304 (1976). 

*2 K. K. Murata and S. Doniach, Phys.Rev.Lett. 29, 285 
(1972). 

''^ N. Buis, J. J. M. Franse, J. van Haarst, J. P. J. Kaandorp, 
and T. Weesing, Phys.Lett. 56, 115 (1976). 
M. Yoshizawa, H. Seki, K. Ikeda, K. Okuno, M. Saito, and 
K. Shigematsu, J.Phys.Soc.Jpn 61, 3313 (1992). 

'^^ M. J. Steiner, F. Beckers, P. G. Niklowitz, and G. G. Lon- 
zarich, PhysicaB 329-333, 1079 (2003). 

*® J. H. Fluitman, Ph.D. thesis. University of Amsterdam 
(1970). 

'^'^ G. G. Lonzarich, Electron (Cambridge University Press, 
Cambridge, 1997), chap. 6. 

L E. Dzyaloshinskii and P. S. Kondratenko, JETP 43, 1036 
(1976). 

'^^ N. R. Bernhoeft, L Cole, G. G. Lonzarich, and G. L. 



11 



Squires, J.Appl.Phys. 53, 8204 (1982). 
^° G. G. Lonzarich, J.Magn.Magn.Mat. 54-57, 612 (1986). 
" A. B. Pippard, Eur.J.Phys. 8, 55 (1987). 
'^^ D. M. Edwards and G. G. Lonzarich, Phil.Mag.B 65, 1185 

(1992). 

J. M. Ziman, Electrons and Phonons (Clarendon Press, 
Oxford, I960). 

The parameters of our model for stoichiometric NisAl at 
ambient pressure are given in Ref. IsrI . It is assumed that 
only the critical parameter a varies appreciably with pres- 
sure p or doping x. The value of a as a function of p 
and X is obtained from its value at p = a; = times 
0(p,x) = Mo(p,a;)/Mo(0,0), where Mo{p,x) is the mag- 
netization in the zero temperature and zero field limit 
(see Equation 0. In the analyses presented in this paper, 
f3{0,x) is inferred from the data in Ref. Is^ and /3(p, 0) is 
taken to be (1 — p/pc) for p < pc (see EquationQand the 
assumed form of a vs p). 

'^^ C. J. Fuller, C. L. Lin, T. Mihalisin, F. Chu, and 
N. Bykovetz, Solid State Commun. 83, 863 (1992). 

^® The data for the sample with the lower Tcurie is only nor- 
malised with respect to Apsf{Tcurie)- Motivated by the 
theoretical simulation Apsf{Tcurie) is determined by a 
T^i'^ fit to the data from 10 K to Tcurie and extrapola- 
tion of this fit to T = K. 

" D. Belitz, T. R. Kirkpatrick, and T. Vojta, Phys.Rev.Lett. 



82, 4707 (1999). 

A. V. Chubukov, C. Pepin, and J. Rech, Phys.Rev.Lett. 
92, 147003 (2004). 

N. Doiron-Leyraud, I. R. Walker, L. Taillefer, M. J. 
Steiner, S. R. Julian, and G. G. Lonzarich, Nature 425, 
595 (2003). 

'° C. Pleiderer, S. R. Julian, and G. G. Lonzarich, Nature 
414, 427 (2001). 

^ C. Pfleiderer, D. Reznik, L. Pintschovius, H. v. Lohneisen, 
M. Garst, and A. Rosch, Nature 427, 227 (2004). 
D. Jaccard, A. T. Holmes, G. Behr, Y. Inada, and 
Y. Onuki, Phys.Lett.A 299, 282 (2002). 
T. Jarlborg, Phys.Lett.A 300, 518 (2002). 

" L. Capogna, A. P. Mackenzie, R. S. Perry, S. A. Grig- 
era, L. M. Galvin, P. Raychaudhuri, A. J. Schofield, 
C. S. Alexander, G. Cao, S. R. Julian, and Y. Maeno, 
Phys.Rev.Lett. 88, 76602 (2002). 

F. RivaduUa, J.-S. Zhou, and J. B. Goodenough, 
Phys.Rev.B 67, 165110 (2003). 

C. Stassis, F. X. Kayser, C.-K. Loong, and D. Arch, 
Phys.Rev.B 24, 3048 (1981). 

G. J. Morgan and J. R. Cooper, J.Phys.F 11, 2091 (1981). 
* N. W. Ashcroft and N. D. Mermin, Solid State Physics 

(Saunders College Publishing, Fort Worth, 1976). 



